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The conductivity of a granular metal or an array of quantum dots usually has the temperature 
dependence associated with variable range hopping within the soft Coulomb gap of density of states. 
This is difficult to explain because neutral dots have a hard charging gap at the Fermi level. We 
show that uncontrolled or intentional doping of the insulator around dots by donors leads to random 
charging of dots and finite bare density of states at the Fermi level. Then Coulomb interactions 
between electrons of distant dots results in the a soft Coulomb gap. We show that in a sparse 
array of dots the bare density of states oscillates as a function of concentration of donors and causes 
periodic changes in the temperature dependence of conductivity. In a dense array of dots the bare 
density of states is totally smeared if there are several donors per dot in the insulator. 



I. INTRODUCTION 

Conduction of samples where metallic granules are sur- 
rounded by some kind of insulator (granular metals) have 
been studied intensively for decades^"®. If volume frac- 
tion of the metal x is large, metallic granules touch each 
other and conductivity is metallic. When x decreases 
and crosses percolation threshold Xc, granules become 
isolated from each other and granular metal goes through 
metal-insulator transition. It is generally observed^"^ 
that at the insulator side of transition the temperature 
dependence of conductivity obeys 



a = ao exp 



1/2- 



Recently, similar temperature dependence was observed 
in doped systems of self-assembled germanium quantum 
dots on the silicon surface^" and CdSe nanocrystal thin 
films^^. Below we talk about both granular metal and 
semiconductor dot structures universally using the word 
dot for brevity. 




FIG. 1. The shape of Coulomb gap in the vicinity of the 
Fermi level. Bare density of states in the absence of long range 
Coulomb interaction is shown by the dashed line. Occupied 
states are shaded. 

In doped semiconductors the temperature dependence 
of Eq. (1) is also widely observed at low temperatures^. 
It is interpreted as the variable range hopping conduc- 
tivity between impurities in the presence of the Coulomb 
gap'''^^ of the density of states (DOS) and is called the 



Efros-Shklovskii (ES) law. In a n-type lightly doped 
compensated semiconductor all acceptors are negatively 
charged and randomly situated, an equal number of 
donors are charged positively. Together all random 
charges create a random potential shifting donor levels 
up and down. This results in finite bare DOS go{p) at 
the Fermi level /i. Long range Coulomb interaction of 
localized electrons then produces the parabolic Coulomb 
gap in DOS 



9{e) = -n'eye^ (2) 



at the Fermi level (Fig. 1) leading to ES law with To = 
(1) C • c^/k^, where C is a constant factor, k is dielectric 



constant and ^ is the localization length. 
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FIG. 2. A sparse array of same size neutral metallic dots 
surrounded by an insulator. 

In contrary to a doped semiconductor, in an array of 
neutral dots the bare DOS at the Fermi level goifJ-) — 
and there is no justification for ES law, which requires a 
non-zero go{p) to begin with. Let us consider a sparse 
array of dots with the same radius R shown on Fig. 2. 
Charging energy levels of such array are shown in Fig. 3. 
Here we assume that dots are large enough, so that spac- 
ing between charging levels / nR is much greater than 
spacing between quantum levels of the dot with a given 
charge. Empty peaks in Fig. 3 correspond to energies 
necessary to charge a neutral dot by the first, second and 
so on electrons transferring them from the macroscopic 
piece of the same metal. Shaded (filled) peaks correspond 
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to, taken with the sign minus, energies necessary to ex- 
tract first, second and so on electrons from the dot, or in 
other words, this is density of states of holes. The Fermi 
level of the array is at zero between two peaks and co- 
incides with the Fermi level of macroscopic piece of the 
same metal. 
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FIG. 3. The bare density of ground states (BDOGS) of a 
clean system of identical dots consists of equidistant peaks. 
Occupied states are shaded. 

We want to emphasize that in each dot we are dealing 
with the ground state at a given number of electrons and 
exclude excited states (higher quantum levels) because 
the ground states of dots determine the low temperature 
hopping transport. Indeed, in the Miller- Abrahams net- 
work^'^ of resistances connecting all dots, the exponen- 
tially large activation factor of each resistance depends 
on probabilities of occupation of a dot by a given num- 
ber of electrons. Exponential temperature dependencies 
of these probabilities, as well as the partition functions 
of dots are determined by ground state energies^'*. Thus, 
the density of states of the array we need can be called 
bare density of gromid states (BDOGS). 

Note that in a lightly doped semiconductor with sev- 
eral equivalent conduction band minima each donor has 
excited states close to the ground state, i. e. situa- 
tion is similar to large dots. Still exponential temper- 
ature dependence of hopping conductivity depends only 
on BDOGS^ 13 

Let us illustrate the role of BDOGS in a clean sparse 
array (Fig. 3). In this case, conductivity requires acti- 
vation of electron-hole pairs or in other words, transfer 
of an electron between two originally neutral dots. Ob- 
viously concentrations of positive and negative dots and 
the hopping conductivity obeys 



a = (TO exp[-Ec/2kBT] 



(3) 



independently of excited states, (Same result can be ob- 
tained in Miller- Abrahams resistor network approach.) 
where = c^/kR is the charging energy. 

How then can we explain observation of ES law? Ap- 
parently BDOGS shown in Fig. 3 shoidd be smeared in 
the vicinity of the Fermi level due to some kind of disor- 
der. 



A simplest source of disorder is distribution of sizes or 
capacitances of dots. Indeed, charging energies of larger 
dots are smaller and this can result in the low energy tail 
of the first empty peak of BDOGS and symmetric high 
energy tail of the first occupied peak. Still in a neutral 
system these tails do not overlap and 50 (m) is zero, so 
that this kind of disorder does not lead to ES law. Sheng 
et al.i assumed that in a reasonably dense system of neu- 
tral granules there is a special distribution of distances 
between granules or their mutual capacitances, which can 
lead to ES law. This assumption was found incompatible 
with other experiments^. But more importantly it was 
noticed® that in a system of large granules made of the 
same metal all granules remain neutral in ground state at 
any distribution of mutual and individual capacitances. 
Therefore, inter-granular excitation of electron-hole pairs 
has a gap. 

This leads to an important conclusion that granular 
metal may have finite BDOGS goin) and show ES con- 
ductivity only if in the ground state of the system gran- 
ules are charged"^"^. Formally one can imagine that this 
happens if granules have different work functions*'^. Sev- 
eral possible mechanisms of such fluctuations were dis- 
cussed in literature. Chui'' suggested that very small 
dots can charge big dots in the case when the former 
are so small that their quantum level spacing exceeds 
the charging energy. We concentrate on the system of 
large enough dots and, therefore, ignore this possibiUty. 
Cuevas et al.^ claimed that even in large dots there are 
large fluctuations of the Fermi level SEp due to random 
positions of neutral impurities in different dots. To our 
mind, this possibility can be rejected using for 6Ep a 
simple estimate of typical fluctuation of the average po- 
tential in the metallic dot due to fluctuation of number 
of impurities there. For three dimensional case one easily 
gets SEp ~ Epc^l"^ I {kpRY/^ , where c is relative concen- 
tration of impurities in the dot, kp is Fermi wave- vector, 
and we assume that size of an impurity is of the order 
of kp^ and its potential U ^ Ep. At large R this en- 
ergy is apparently smaller than charging energy / kR 
even at c ~ 1, so that dots remain neutral. Baskin and 
Entin^ considered fluctuations of surface part of the dot 
energy as a reason of ionization of dots. Again if we as- 
sume that these fluctuations are result of random distri- 
bution of point-like impurities located on the dot surface 
we come to conclusion that corresponding fluctuation of 
energy decreases like 1 /R^ and can not compete with the 
charging energy. Thus, simplest internal mechanisms of 
fluctuations of the work function mentioned above are 
too weak to charge array of three-dimensional dots. 

In this paper we study models of arrays of dots af- 
fected by external disorder, where the origin of charging 
and BDOGS go is transparent and in some cases con- 
trollable. We assume that the insulator between clean 
dots has a concentration of donors No with electron en- 
ergy Epi higher than Fermi energy /i in dots, so that at 
low temperatures all donor donate an electron to dots 
and charge them. For a given concentration of dots N 
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(Fig. 2) wc can introduce the average number of extra 
electrons per dot v = N^/N, which by analogy with 
quantum Hall effect can be called the filling factor. Max- 
imum number of electrons which can be added to a dot 
is Umax — {Ed — jj) / {e^ / kR) . We assume that u < Umax 
so that all donors lose their electrons. In the semicon- 
ductor language we are dealing with compensated p-type 
semiconductor where dots play the role of multi-charged 
acceptors. 

In Sec. II we summarize our main results for BDOGS 
at the Fermi level gaip) as function of the filling factor 
I/. In Sec. Ill the sparse 3D array of dots (Fig. 2) is 
discussed. In Sec. IV we extend this discussion to 2D 
array. In Sec. V we study in detail the super-dense 3D 
array shown in Fig. 4, which is the other extreme, and 
then comment on the realistic moderately dense array 
(Fig. 14). We also see how the properties of sparse 3D 
arrays cross over to dense arrays. After studying differ- 
ent arrays one can see that external doping brings about 
non-zero 5o(m) and leads to ES law. 
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FIG. 4. A "super-dense" array of dots. Only one donor is 
shown by + and electron donated by it is shown by -. 
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FIG. 5. BDOGS g^is) at a certain filling factor. Occupied 
states are shaded. 

As a result 50 (a*) may oscillate with v. For example, at 
1/ = 1/2 the Fermi level \x is in the middle of the BDOGS 
peak and for that reason gaip) can be rather large. On 
the other hand, at = 1 the Fermi level /i is in a tail 
between two BDOGS peaks, and 50 (m) tends to be much 
lower. Actually the total dependence of 50 (/x) on f is 
somewhat more complicated and consists of three parts, 
the growing part, the oscillating part and the saturation 
part (Fig. 6). The first (growing) part takes place at 
small v, where the situation is similar to p-type semicon- 
ductors at low degree of compensation^^. We get BDOGS 
growing linearly with v 



5<(/x)~«i/iVV3/e2 



(z/ « 1). 



(4) 



In the second (oscillating) part, fi typically dwells in 
BDOGS peaks, 



(5) 



II. SUMMARY OF RESULTS 



We start from sparse three dimensional arrays of dots 
of the same radius R, which are randomly situated in 
space with concentration N <C (47ri?3/3)~^. We assume 
that the dots are big enough so that Coulomb effect over- 
weights quantum level spacing. For sparse arrays of dots, 
doping introduces two types of charges: positive, empty 
donors and negatively charged dots. Both are randomly 
situated and create random potentials growing with v. 
These charges result in two effects (Fig. 5). Firstly, 
the (5-shaped peaks in BDOGS go{e) become somewhat 
smeared, since the energy it takes to bring an electron 
to a dot is affected by the random potential (the effect is 
similar to that of gedanken random gate potential). As 
a result, each peak gets tails. Secondly, electrons com- 
ing from donors fill some dot states and hence move the 
position of the Fermi level up. 



When V is very close to integers, the Fermi level drop 
into minima of go{£) where 

j,8/3 JV5/3^2 

Qminill) fa/(^^^ii) (1< < ys). (6) 

Here Vs ~ l/(iVi?3)i/4 >> 1 is the filling factor at which 
oscillations become relatively small and the third (satu- 
ration) part starts over. This happens because fluctua- 
tions of the Coulomb potential are so big that BDOGS 
is almost uniform everywhere 



The number of large oscillations is big only when dots 
are far from each other. Here and below we often omit 
numerical coefficients. 
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FIG. 6. BDOGS at the Fermi level (?o(a*) of a sparse 3D ar- 
ray as a function of filling factor (solid line). The reference to 
the equation appropriate for a part of the curve is shown next 
to it. Dashed lines describe locations of minima and maxima 

of the oscillating part. 

For less ideal arrays where the dots have shghtly dif- 
ferent sizes, the distribution of sizes can wash away the 
oscillations to a certain extent, making the line of max- 
ima lower and the line of minima higher. As a result the 
two lines approach each other faster {i/g is smaller) than 
the situation of identical sizes. 

Oscillations of 50 (a*) lead to periodic transitions be- 
tween the ES law and the Mott's law at a given low tem- 
perature as shown in Fig. 7. This happens because in the 
very vicinity of fj, the long range Coulomb interaction cre- 
ates the parabolic Coulomb gap (not shown in Fig. 5). 
The width of Coulomb gap as follows from Eq. (2) and 
Fig. 1 depends on go 



(8) 



At large go the Coulomb gap is wide (A >> ksT) and 
the conductivity obeys the ES law (Eq. (1)), which does 
not depend on 90- At a very small go the width of the 
Coulomb gap is smaller than ksT and becomes irrelevant 
for conductivity. In this case we arrive at the Mott's law 
region 



a = ao ■ exp 



/5 



go{li)a^T 



1/4- 



(9) 



with strong dependence on goifj)- Here /? is a numerical 
coefficient'' . 

Oscillations of BDOGS 50 (m) with u lead to oscillations 

of other measurable quantities. It was shown "'^ that in 
systems with the Coulomb gap BDOGS 50 (a*) plays the 
role of the thermodynamic density of states dn/dfi, where 
n is concentration of electrons, so that it determines 
screening radius of the system — {ATrgoe^ / n)^^^^ used 
below for self-consistent calculations of BDOGS. One can 
also measure BDOGS at the Fermi level with help of ex- 
tremely low temperature (lower than typical quantum 
level spacing) specific heat and microwave absorption^' 
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FIG. 7. Ranges of ES and Mott's laws alternating with 
growth of filling factor v at a given low temperature. 

Until now we talked about 3D arrays of dots. In Sec. 

IV wc consider arrays of dots located in a 2D plane. In a 
2D array dots can be charged by donors located in paral- 
lel to the plane (^-layer)^*^'^^. In 2D case, however, there 
is a more practical way to charge dots using a metallic 
gate parallel to the plane of dots. At >> 1 results 
for (7o(m) are almost independent on the way of charg- 
ing, because most of the random potential is created by 
dot charges. Dependence of 50 (m) on v in 2D qualita- 
tively looks similar to the 3D case discussed above, but 
quantitatively it is somewhat different (see Sec. IV). 

Dependence of 50 (a*) on v in 2D should leads to oscilla- 
tions of conductivity. Some oscillation were observed in 
experiments^°'^^. But there are additional ways to mea- 
sure this dependence in a gated in 2D structure. First^'^, 
one can study change of conductivity where the gate volt- 
age charges dots. Second^®, one can study inverse small 
signal AC capacitance of unit area 1/C, which is related 
to the two-dimensional screening radius of the dot system 
Ts = K/{2Trgoe'^) by equation 



1/C = 4TrK{do+rs), 



(10) 



where do is distance from gate to the plane of dots. This 
method is similar to magneto-capacitance measurements 
in quantum Hall effect^". 

Until now we talked about sparse arrays of dots. For a 
sparse array of dots the variable range hopping conduc- 
tivity is measurable only when the tunneling length in 
insulator is large enough. This probably can be achieved 
in some cases in narrow gap semiconductors. Otherwise 
conductivity is so small that it can not be measured. 
In this case, thermodynamic measurements of oscillating 
quantities may be more convenient. 

It is much easier to measure conductivity in a dense 
array of dots. Actually most of cited experiments are 
done at NR^ ^1. In Sec. V we study clean metallic 
cubes separated by thin layers of an insulator doped by 
donors in the "super-dense" limit d « R (Fig. 4). It 
is assumed that the tunneling conductance between two 
cubes G << e^/h. We show that in the dense array 
BDOGS as a function of the average number of donors 
per dot u behaves as shown in Fig. 8. At << 1 BDOGS 
is very small and grows super- linearly with u (like v^), 
and at 1/ ^ 1 it saturates at the value gs ^ K/{e'^RDd). 
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FIG. 8. BDOGS at the Fermi level godj-) of a "super-dense' 
array as a function of filling factor v. 



III. BDOGS IN SPARSE THREE-DIMENSIONAL 
ARRAYS. 



1. Saturation of BDOGS 

With enough impurities the BDOGS is strongly 
smeared and it approaches an averaged value given by 
Eq. (7) . This happens when the fluctuation of potential 
7 is large enough, namely 

7 > e^/KR . (11) 

In order to find the critical i^s of the saturation 
BDOGS, let us study the fluctuation of electric potential. 
For a given quantum dot, all the charges within a dis- 
tance of the screening radius r., contribute into a collec- 
tive Coulomb potential. There are two kinds of charges, 
positive donors and negative dots. Let us compare the 
charge fluctuations made by them. The average number 
of impurities within is Nur^, with a typical charge 
fluctuation e^^Nur^. For the charged dots, the typical 
fluctuation of dot number in is \/iVrf , and each dot 
carries a charge o f ei^ o n average, so the typical charge 
fluctuation is ve^jNrl. At >> 1 the charge fluctua- 
tions produced by charged dots are larger than that of 
donors, so that the fluctuation of the potential energy 
of an electron is mostly related to charge fluctuations of 
dots 




On the other hand, the BDOGS at the Fermi level 
5o(a*) determines the linear screening radius due to re- 
distribution of electrons between dots 

1 1 

rs = , !=■ (13) 

v/4^i^C«)eV« 

With the help of Eqs. (12) and (13), the require- 
ment (11) becomes 

7 ~ e^v{N/Rf'^/K > b'^/kR . (14) 



Therefore only at v > Vs = 1/{NR'^Y/'^ can we have 
big enough potential fluctuation. M v > Vs BDOGS is 
smeared so strongly that at a given v BDOGS g{e) ~ Qs 
at all energies, and 50 (a*) does not depend on the position 
of n any longer. 

2. Line of maxima of oscillating BDOGS. 

When V is close to a half-integer M — 1/2 (M is an 
integer) the Fermi level is in the middle of the Mth peak 
in the BDOGS of dots (Fig. 9). 




FIG. 9. The origin of Qmax- The solid horizontal line is 
the Fermi level, the dashed line is the fluctuating potential 
energy. Two peaks of BDOGS of the array are shown, as well 
as two corresponding levels of 4 dots. 



The single dot charging levels vary with the electric 
potential energy. The typical fluctuation of the potential 
energy is determined by Eqs. (12), which broadens the 
BDOGS peak (Fig. 9). The BDOGS at the peak is 

N 

Qmaxil^) ^ 9peak{s) • (15) 

7 

In such a situation the linear screening radius is 



max 

At 1 < V < Vs Eqs. (12), (15) and (16) self-consistently 
determine all three unknowns: the BDOGS Qmax^lj)-, the 
screening radius and the typical fluctuation of the poten- 
tial energy. We arrive at Eq. (5) and 

rs ~ !/2/3/^l/3 » ^-1/3 ^ (17) 

7~eV/37VV3/^. (18) 

It is clear from Eq. (5) that Qmax arrives at the value 

Qs ^ N/{e^/KR) when z/ - i/,. = 1/iNR^y/* » I. Note 
that the theory of this subsection is similar to Ref.^-^. 

3. Line of minima of oscillating BDOGS. 
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FIG. 10. The origin of Qmin- The meandering hne rep- 
resents the fluctuating potential energy as a function of co- 
ordinate. The regions of the shortest size ~ I where dots 
are completely filled with electrons or holes arc shown with 
thicker lines. The solid horizontal line in the middle is the 
Fermi level, the other two horizontal lines correspond to the 
nearest peaks of an isolated dot. Two peaks of BDOGS of 
the array are shown, as well as two corresponding levels of 7 
dots. 

When v is close to an integer M the Fermi level lies 
in the tail between the Mth and (M -I- l)th peaks in the 
BDOGS of dots. Due to the low BDOGS at tails, screen- 
ing of the long-range fluctuations is weak and strongly 
nonlinear. The random potential fluctuation grows un- 
til it reaches /2k,R, where levels of some dots cross the 
Fermi level (Fig. 10). This picture is similar to the model 
of compensated semiconductor and we treat it below fol- 
lowing ideas of Chapter 13 of Ref.^. Let us define rs 
as nonlinear screening radius. Then all space can be di- 
vided into blocks of size with fluctuating excess charges 
q r-j ue^jNr'j which create potential 7(rs) ~ ^^/Nr^. 
This requires 



eVKi?~7(rs) ^/Nrf (19) 

and we arrive at 

The potential fluctuation forms wells and humps 
(Fig. 10). Electrons fall into wells and holes fall into 
humps. We will discuss the wells and electrons, although 
our discussion is equally applicable to humps and holes. 
Inside a well of size rs there are other humps and wells 
of shorter range. Electrons inside wells find themselves 



in wells of smaller scale and so on until they reach the 
shortest scale I. 

The shortest range / can be cletermincxl by the fact that 
each dot accepts only one electron. The exce ss nu mber 
of charges inside such a well is of the order vV Nl^. The 
maximum number of electrons in the well is of the same 
order, otherwise the well turns into a hump. On the other 
hand the maximum number of electrons is Nl^ so that 
Nl^ ~ i/y/NP. This gives 

I ~ {v'^/Nf/^ (21) 

for the characteristic size of the shortest scale well, and 

^il) ^ ^^TVP ~ Z/4/3^1/3g2/^ ^ (22) 

for the characteristic depth. It can be checked that 
I « Ts as long as << Vg, in accordance with our 
assumption. 

As we said, in a cube of si ze Ts, the fluctuation of ex- 
tra charge is g ~ ve^J Nrf, which is balanc;c;(l by the 
screening charge ANr^e {AN is concentration of dots 
participating in screening): 

AiVr^e ~ uey/Nrl , (23) 

therefore 

AN - v'^N'^R^ . (24) 

With the help of Eqs. (22) and (24), one can estimate 
the BDOGS at the Fermi level as gmmi^^) ^ AN/j{l), 
and arrive at Eq. (6). For peaks oi g{e) at e = ±e^/2KR 
we have g{e) « gs » grninip) (see Fig. 10). 

At ly ^ Vs the shortest scale of wells and humps / ap- 
proaches the longest scale of wells and humps rs, and 
each dot can accept more than one electrons or holes 
because 7 > c^/kR. Therefore the non- linear screening 
crosses over to linear screening at v ^ i>s- At the same 
time the fluctuation of potential energy 7 is big enough 
to smear the BDOGS. It can be checked that gmin and 
grnax mcrgc to gs &t V ^Vs. 

4. Grov^fing part of BDOGS. 

When concentration of impurities is rather low, v « 
1, an electron released by a donor would like to stay on a 
dot close to the donor. The common situation is a donor 
accompanied by one charged dot (1-complex), while some 
donors stay alone (0-complex), and some donors are ac- 
companied by two charged dots (2-complex). Charge 
conservation requires A^o = A2, where Nq and N2 are 
concentrations of 0- and 2-complexes. 

At !^ << 1 nearly all donors are independent, and A^o, 
Ni and N2 are proportional to v. Therefore the Fermi 
level remains a constant and the BDOGS grows linearly 
with V sX u «1. 

The details of this problem is the same as in a weakly 
compensated doped semiconductors. It has been stud- 
ied quantitatively^^, counting the numbers of 0- and 2- 
complexes. The quantitative results are 
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Nq = Nv exp 



N2 = 7.14 • 10-4 • ( ^ 



V 3 



■Nv 



2_/yl/3 
KjJL 



(25) 



(26) 



(27) 



As the concentration of impurities Nv — Nq + Ni + N2 
is not changed, the concentration of electrons on dots is 
controlled by the Fermi level /z, therefore 

^dN2_dNo^ 0.2nvN^/ye^ . (28) 

We kept all coefficients in this subsection because they 
are known. 

Thus g<:{lJ-) is proportional to v and follows Eq. (4). 
This picture works well a.t v « 1. At v > 1 each dot 
is typically affected by more than one impurities. It is 
straightforward to check that Eq. (4) matches Eq. (5) at 
V r^l. 

5. Periodic array of dots. 

Above we have discussed sparse arrays in which the 
dots are situated randomly. Now we turn to a question, 
what if the dots arc arranged periodically in space? 

The qualitative picture in Fig. 6 is still mostly correct. 
At >> 1 there are still the oscillations of 50 (a*) with 
V. One should notice that the strongly charged periodic 
dots do not contribute to the fluctuations of electric po- 
tential and only donors produce the fluctuations. As we 
discussed above these fluctuations are weaker than in the 
random sparse array. Therefore more donors are needed 
to smear BDOGS to the same extent. In other words, the 
lines of maxima and minima in this case converge slower 
and the value of Vg is larger. 

The linear growth of BDOGS at << 1 shown in 
Fig. 6 is not valid in a periodic array, since 0- and 2- 
complexes originate from the random distribution of dots 
in space. In periodic array a typical donor donates an 
electron to the nearest dots and forms an electric dipole. 
It takes additional energy to take an electron out of a 
dipole and put it near another dipolc. Therefore a hard 
gap is formed at the vicinity of the Fermi level. Only the 
rare cases where several donors arc close to each other 
can produce non-zero BDOGS at the Fermi level. As 
a result goifJ') does not grow linearly with v, but with 
higher power. The similar situation will be studied more 
thoroughly in the dense array of Sec. V. Therefore, here 
we gave only a brief description of the result. 



IV. BDOGS IN A SPARSE TWO-DIMENSIONAL 
ARRAY 



The 5o(a*) as a function of v in 2D looks similar to 
3D (Fig. 6). It also has growing part, oscillating part 



and saturation part. However, in 2D oscillations survive 
longer than in 3D, because the collective potential in 2D 
is much weaker. Indeed, the fluctuation of potential en- 
ergy in 3D grows with Vg as a power law 



73D 



-^jNr^ DC r 



1/2 



(29) 



while in 2D the long-range fluctuation of potential energy 
grows slower than logarithmically with 

72Z? ~ j C^fNrdr oc ^JM}/^)- (30) 



JV 



Therefore the line of minima of oscillating part in 2D 

is quite different from that of 3D. The tail in the g{e) vs. 
e curve is produced here by Coulomb interaction of near- 
est neighbor dots, and the minimu m value of g appears 
to be (at KKu <Vs = l/\/NR?) 

9min ~ N^{ve''/Kf/e |e^eV(«iJ) ~ v'^N^R^k/c^. (31) 
The 2D linear screening radius produced by Qmin is 
K 1 



2'ng„ 



(32) 



One can verify using Eq. (30) that 72D << / kR for 
such and therefore long-range contribution can be ne- 
glected. 

The growing part at < 1/2 can be studied in the 
similar way as in 3D. That is, one can use the idea about 
0-complex, 1-complex and 2-complex to find the Fermi 
level 

II ~ -e^N^/'^/n (33) 
and BDOGS at the Fermi level {aX u « 1) 

5<(/x) ~ KvN^'^/e\ (34) 

For the maximum line of oscillating part, the collec- 
tive typical fluctuation of potential energy is given by 
Eq. (30) and BDOGS at the Fermi level is 

TV N'^/'^K, 

9max{^J) ~ — l«V<Vs = 1/VNR'^. 
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(35) 



Here we dropped the logarithmic factor and therefore 
disappeared from the equations. In this approximation 

unlike 3D situation, we do not need the self-consistent 
method. As a result gmaxil^) decreases with v much 
slower than in 3D. 
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Both Eqs. (35^ind (31) reach .g,(M) ^ N/i^/nR) at h = WPio QiQj = J + J E^^^ Q^^^' • 



1/ ~ = 1/\JNR?, where the distinction between peaks 
and valleys disappears. 

Up to now we discussed charging of two-dimensional 
array of dots by donors or a gate situated outside of dots. 
In some two-dimensional systems electrons may be pro- 
vided to a dot by donors located close to the dots. This 
can be done, for example, by creating a two-dimensional 
gas with help of close layer of donors parallel to its plane 
and then by etching outside dots both the gas and the 
donor layer. If such dots are essentially two-dimensional 
and heavily doped, internally induced fluctuations can 
cause their charging, in contrary to what we said in in- 
troduction about three-dimensional dots. Indeed, in this 
case, fluctuations of number of electrons in the dot and 
correspondingly of the Fermi level decrease inverse pro- 
portionally to the dot radius. This decay is similar to 
the charging energy, so that ratio of these energies is just 
proportional to the ratio of kinetic and potential energies 
of electron in the dot^^. 



V. DENSE THREE-DIMENSIONAL ARRAYS OF 
DOTS 

Until now we dealt with sparse arrays of dots where 
NE? << 1 as shown in Fig. 2. Let us discuss what 
happens when the small parameter NE? grows. The sat- 
uration filling factor = (NR^)^^''^ decreases, which 
means the number of oscillations decreases as well. The 
oscillations disappear at i^g ^ 1 when NR^ 1, i. e. 
when the typical separation between dots is of the order 
of the dot size. 

Let us now concentrate on the other extreme (Fig. 4) 
where say cubic dots are closely packed in a cubic lattice 
with the period R and the width d of an insulator be- 
tween dots is small {d « R). The insulator is uniformly 
doped by donors. We show below that at a large enough 
average number of donors per dot, v > 1, the BDOGS is 
smeared. On the other hand at << 1 typical donors 
do not contribute to BDOGS at the Fermi level but rare 
clusters of donors result in a small BDOGS. 

For the system of metallic cubic dots {i} separated by 
a clean insulator, potentials (pi and charges Qi of dots 
are related by linear equations 



and 



where Pi 



-Py Qj 



(36) 



(37) 



ij is the reciprocal matrix of the capacitance ma- 



trix Cij . The coefficients satisfy Pij = Pji , < Pij 
and depend on the relative positions of i and j only. 
Without donors the Hamiltonian can be written as 



(38) 



Of course, it has minimum when all Qi = 0. 

In order to understand the role of ionized donors let 
us study an isolated donor between the dots A and B 
(Fig. 11). The positive donor is completely screened by 
negative charges qa and qb appearing of the surface of 
A and B, i. c. qa + Qb + e = 0. These charges leave 
compensating positive charges —qa and — in the dots 
A and B. If a;^ and xb are distances from the donor to 
the dots A and B respectively then 



QAXB = QbXA- 



(39) 



One can easily arrive at Eq. (39) imagining that the donor 
for a moment is replaced by the uniformly charged plane 
parallel to the faces of cube A and B (with the same xa 
and xb), and minimizing the sum of Coulomb energies 
of two emerging plane capacitors with distances between 
planes xa and a:;^ as a function of qa- Note that at 
d « i? all points of the charged plane are in the same 
conditions and the plane fields arc just a superposition of 
fields of these charges. This means Eq. (39) should also 
hold for a point charge of this plane or just for a donor. 
Eq. (39) means that potentials of dot A and B are equal. 
In other words, the group made of the positive donor, 
qa and qb produces an electric field which is confined in 
the vicinity of the donor. The group does not affect the 
electric potential of any dot. Therefore the group can 
be totally forgotten when we study the charges and po- 
tentials of the dots, or in other words it can be totally 
ignored in the Hamiltonian. 




< 


> 







FIG. 11. Two-dimensional cross section through a donor 
(with charge +e) located between the dots A and B for the 
system shown in Fig. 4. The group within the dashed circle 
is neutral altogether and docs not affect the energy levels of 
the dots A and B. The dot A gets an electron donated by the 
donor because the donor is closer to its border. 
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On the other hand, the compensating positive charges 
—QA and —qb can not be ignored. One can say that the 
donor's charge is split into two fractional parts —qa and 
—qb- Charges —Qi are uniformly distributed between 
and e. 

When there are more than one donor around a dot A 
we can simply add them up Qad = —QAi — 5A2 — • • • to 
get the total positive charge induced by donors on the 
dot A. 

Each donor also donates an electron to the array of 
dots. We assume that d » a, where a is the tunneling 
length of electrons in the insulator so that the conduc- 
tance G between two dots is much smaller than e'^/h. In 
this case each dot can get only integer number of extra 
electrons Ui . The neutrality of the whole system requires 
— e^rij = '^QiD- For a given fractional set {Qio} the 
integer set {n,} minimizes the Hamiltonian 



is the effective dielectric constant of the array. For 
charged dots stability of the global ground state requires 
e^^) > /z = and e\°^ < /x = 0. 

Let us study the BDOGS in the case when there are 
more than one donors per dot (i^ >> 1). In order to 
find BDOGS we neglect the interaction terms in Eq. (40) 
and consider only diagonal terms. As a result the ground 
state of the system is realized by such a set of {iii} that 
—e/2 < QiD — me < e/2. So effective charges on differ- 
ent dots in the ground state of the array are uniformly 
distributed from —e/2 to e/2. This gives the energy in- 
dependent BDOGS 



H 



riic) + 



2 I \ 



^ 2 



{Qw - nie){QjD - rije) 



(40) 



This energy is calculated from the reference point of the 
sum of energies of neutral groups around each donor, 
which does not depend on n,. If each donated electrons 
were shared by dots this fractional charges would neutral- 
ize QiD and make H = and all cpi = 0. This is equiv- 
alent to what would happen if all grain are grounded. 
Discreteness of rii makes finding the ground state rii non- 
trivial. It is clear that in the globally neutral system the 
chemical potential jj, lies in the same place as for a system 
of neutral dots, i. e. /u = 0. In the ground state we can 
introduce "one-dot" energies of an empty dot state above 
the Fermi level assuming that charges of other dots are 
fixed 



1 



^ = H{n, + 1) - Him) = -ecPi + ^Pue^ 



where 



(41) 



(42) 



is the electrostatic potential of the dot i. For an occupied 
state in the similar way 



4°^ = Him) - Him - 1) = -e<l>i - -Pue^ . (43) 



In the absence of donors when all dots arc neutral 



1. e. 

0, we get s\'" = +^Pue' and ef> = -\Pue\ The 
single electron BDOGS for neutral dots looks the same 
as in Fig. 3 except the separation of peaks or the level 
spacing is smaller now e^Pu <^ e^/KeffR — (the 
calculation of Pu is in Appendix). Here 



9s{e) 



1 



e-'PiiB? 



(45) 



Consider now the role of non-diagonal terms. At a big 
distance » R the coefficients Pij have the Coulomb 
form Pij « 1 /Kef fVij . They, therefore, result in a density 
'of ground states (DOGS) with the Coulomb gap at the 
Fermi level, which in turn leads to ES law. 

The finite BDOGS at the Fermi level and ES law at 
v » 1, originates from random distribution of charges 
Qi, which in turn is related to the fact that several donors 
contribute into each Qio- Even for nearest neighbor dots 
QiD and Qjo are independent since they include contri- 
butions of different donors. 

We turn now to the case of a small density of donors 
y « 1 and first deal with typical donors isolated from 
each other. As discussed above for such an isolated donor 
between the dots A and B, the charges induced by the 
donor are related by Qad + Qbd = +e. If the donor is 
closer to the dot A, Qad > e/2 > Qbd > 0. If wc ne- 
glect interaction between the dots A and B in Eqs. (41) 
and (43) we arrive at finite BDOGS of isolated donors 
at the Fermi level. States at the Fermi level are pro- 
vided by dots where Qbd = 1/2 — S with 5 « 1. Then 
the dot B has a state just above the Fermi level and 
A has a state just below it. The fact that these states 
are so closely correlated in space leads to their elimina- 
tion by the interaction between the dots A and B. In- 
deed, in the groimd state the electron donated by the 
donor is on the dot A. The total charges on A and B is 
Qa = Qad - e = -Qbd and Qb = Qbd- The electric 
potential on B is (j)B = PbbQb + PbaQa and according 
to Eq. (41) all the levels on B are moved down by e(j)B- 
Using the fact < Qbd < e/2, we have 



< e</.B < (1 - a)PBB eV2 , 



(46) 



nR/d 



(44) 



where a = Pba/Paa = 0.34 ± 0.01 represents the in- 
teraction between nearest neighbor dots (see the calcula- 
tion of a in Appendix). According to Eq. (41) the low- 
est vacant level of dot B is at least aPae^ /2 above the 
Fermi level. Similarly, according to Eq. (43) the high- 
est filled level on the dot A is moved up by —e(j)A with 
< —e4)A < (1 — Q;)PAAe^/2 and it is at least aPae^ /2 
below the Fermi level (Fig. 12). 
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FIG. 12. BDOGS of isolated donors at v « 1. Energy e 
is measured in units of aPae^ /2. There is a hard gap of the 
width aPiie^ around the Fermi level. 



Thus, wo showed that for z/ << 1 typical isolated 
donors have a substantial hard gap of width aPue^ in 
their density of states. This gap can not be destroyed by 
long range interactions of such donors because they form 
neutral complexes and weakly interact with each other 
at long distances. 

One can also think of a donor in the neutral complex 
with electron of the neighboring dot A as a hydrogen-like 
donor in a semiconductor. Then the hard gap obtained 
above is similar to Hubbard gap. The electric dipole 
made of a donor and an electron in principle can accept 
another electron with a weak binding energy, but the 
energy difference or Hubbard gap between the two elec- 
tronic levels makes sure that only the lower level is filled 
in the ground state. 

Up to now we talked about typical isolated donors at 
V « 1. In fact there are rare configurations that make 
a non-zero BDOGS at the Fermi level. Let us show how 
this happens. 



1/5 



1/5 



1/5 



1/5 



1 r 



1 



FIG. 13. An example of rare configurations leading to the 
finite BDOGS at the Fermi level at v « 1. The cluster of 
5 dots and 4 donors is effectively positively charged in the 
ground state. Charges of 5 dots are shown in units of e. All 
the other dots in the vicinity of the cluster remain neutral. 

Consider, for example, a cluster of 4 donors around 
one dot (Fig. 13). The charges of donors are shared by 



5 dots in such a way that charge of each donor is effec- 
tively split into 4e/5 going to the central dot and e/5 
going to its nearest neighbor. The middle dot has effec- 
tive charge QiD = 16e/5 and each of the four neighbor 
dots has QjD = e/5. Let us show that in the ground state 
of the cluster = 3 for the central dot and rij = for all 
others, i. e. there are only three electrons on the cluster 
of four donors so that the cluster has a net charge +e. 
It is effectively shared by the 5 dots, each getting +e/5. 
According to Eq. (41) addition of the fourth electron to 
this cluster costs 



-(e) 



-1(4.- 



1) 



P 



(47) 



where a = P12/P11, and as shown in Appendix a = 

(e) 

0.3405. This means that el ' > ^ = and the positively 
charged cluster of Fig. 13 is stable with respect of bring- 
ing the fourth neutralizing electron. The physical reason 
for such stability is uniform smearing of the net charge 
over 5 dots which diminishes potential of the cluster (f>i 
at the central dot attracting the fourth electron (the first 
term in the right side of Eq. (47)). 

We have shown above an example of a positive "donor- 
like" cluster. Moving each of the four donors away from 
the central dot in the direction of the neighbor dot one 
can construct an "acceptor-like" cluster in which charge 
— e is evenly shared by 5 dots. Such a cluster accepts an 
electron released by the cluster in Fig. 13. Thus, some 
clusters charge themselves by "self-compensation" sim- 
ilarly to amphoteric impurities in semiconductors. Of 
course, there are also neutral clusters and tuning posi- 
tions of donors between dots one can continuously go 
from a neutral cluster to the charged one. This means 
that BDOGS at the Fermi level created by clusters is 
finite. 

Using coefficients Pij calculated in Appendix, it is easy 
to find that the smallest charged cluster in the ground 
state is the cluster of 4 dots with 3 donors between them. 
Therefore, at v « 1 the BDOGS at the Fermi level 
grows as 50 (m) oc and it begins to be significant at 
1/ ~ 1 (see Fig. 8). 

The structure of self charging clusters discussed above 
is similar to the mechanism which provides a non- 
zero bare density of states in a system of quasi-one- 
dimensional electron crystals at low impurity concentra- 
tion-^"^. There is also similarity to self-compensation of 
clusters of several donors in uncompensated semiconduc- 
tor near the insulator-metal transition^*. 

As we saw above in the "super-dense" array of dots 
with << i?, a large enough concentration of donors 
can easily smear BDOGS. This is a natural explana- 
tion of the experiments, where the ES law is observed. 
Let us consider parameter Tq of the ES law (Eq. (1)) 
for a dense array of dots at v » 1. According to 
Ref.* To ~ e^/{Keff^), where ^ is the localization length 
for tunneling to distances much larger than R. When 
an electron tunnels through insulator it accumulates di- 
mensionless action d/a, where a is the tunneling decay 
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length in the insulator. On the distance x electron ac- 
cumulates x/R such actions. Thus, its wave function 
decays as (2xp{—xd/ Ra) = cxp(— x/^), where the local- 
ization length ^ = aR/d. This enhancement of localiza- 
tion length is similar to the one derived for disordered 
granular system near percolation threshold^^. It leads to 



KeffttR KR'^a 

We see that with decreasing d the temperature To de- 
creases as d"^. This continues until conductance of the in- 
sulator layer G — (e^/?i)(i?A;F)^e~^''/° reaches the quan- 
tum hmit G ~ e^/?i, i. e. while d > dc = a\n{RkF)- At 
d « dc the characteristic ES temperature To(d) should 
vanish when d ^ 0, but the way how this happens is 
still unknown (see also recent publications^^'^'' concerned 
with the case G » 1). 

In the above calculation of the localization length we 
assumed that tunneling through a metallic dot does not 
accumulate more action. Actually a dirty dot provides 
logarithmic contribution to the tunneling action^*. We 
assume that it is much smaller than 2d /a resulting from 
crossing the insulator. 

Until now we discussed only the "super-dense" array 
shown in Fig. 4. Let us consider what happens in the sys- 
tem of almost densely packed metallic spheres (Fig. 14), 
where insulator occupies a larger fraction of space than 
in Fig. 4. 



d 




FIG. 14. An almost densely packed array of metallic 
spheres {d « R). A donor is shown by + and electron 
donated by it is shown by -. 

One can show that in such a system BDOGS as a func- 
tion of number of donors per dot v behaves qualitatively 
similar to Fig. 8. At smaller i> only a cluster of donors 
around a dot can create a finite BDOGS, while at ly >> 1 
donors easily smear BDOGS. In this case, characteristic 
temperature of ES law is Tq w e^d/aR. It loses one 
power of d/R because the dielectric constant of such sys- 
tems does not diverge at d ^ 0. When distance between 
dots d grows and exceeds R, we get sparse periodic ar- 
ray and recover oscillations of BDOGS at >> 1. At 
u << 1 there is a smooth cross-over between two super- 
linear growths of BDOGS related to rare donor clusters 
in the vicinity of some dots. 

Although we arrived at >> 1 as the universal cri- 
terion of substantial smearing of BDOGS for both dense 



arrays (Figs. 4 and 14) and sparse arrays(Figs. 2), we 
should understand that the required concentration of 
donors Nd in the insulator grows while volume fraction 
of insulator decreases from the arrays on Figs. 2 and to 
the one on Fig. 14 and then to the one on Fig. 4. 

In some cases metallic granules are first coated by the 
layer of a doped insulator (for example, the metal's oxide) 
and then compressed into the bulk array with a clean in- 
sulator of different kind of filling or just air filling empty 
space (Fig. 15). For simplicity we assume the dielectric 
constant of the coating insulator is close to that of the 
filling space between dots. 




FIG. 15. An array of coated dots. The width of coating 
insulator is ui. A donor is shown by + and electron donated 

by it is shown by -. 

We would like to consider BDOGS for such a model 
of granular metal in order to understand cross-over be- 
tween the studied above case of donors residing outside 
dots and those on the dot surface or inside dots. In the 
latter case, donors are screened and act as short range 
neutral impurities. According to the introduction such 
impurities can not charge a neutral array of dots. Thus, 
effect of donors in the layer of the width w « R should 
decrease as w vanishes. 

Indeed, one can calculate, fluctuation of average the 
potential of donor layer of the width w which plays the 
role of fluctuation of work function of dots. It can be es- 
timated as a potential created by the layer of the width w 
by fluctuating concentration [NjjwR^Y^^/wFi^. One can 
imagine that, correction to the work function is created 
by effective random double layer. Such correction to the 
work fimction is of order {{Ndw^Y^^ / kR. It becomes 
larger than the charging energy only at Ndw'^ » 1^^. 
This condition of BDOGS smearing is obviously much 
stronger than "one donor per dot" condition valid for a 
uniform insulator. When w approaches the lattice con- 
stant the role of donors is so much weakened that in order 
to smear BDOGS one needs the relative concentration of 
donors to be of the order of 50%. At the surface, donors 
and their images arc still forming dipoles with fluctuat- 
ing dipole moments, which create fluctuating double layer 
potential. On the other hand, when donors enter inside 
the dot, each donor becomes exponentially screened from 
all sides, the dipoles and the random double layer disap- 
pear. This diminishes the role of donors further so that 
as we said in introduction even large concentration of im- 
purities inside the whole three dimensional dot can not 
charge the dot. 
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VI. CONCLUSION 
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This paper addresses the problem of explanation of the 
temperature dependence of the conductivity observed in 
granular metals and arrays of quantum dots. We show 
that ES law is the most natural explanation for this de- 
pendence (in spite of recently expressed doubts^"). For 
this purpose wo present a simple model with random 
charging of clean metallic dots in both sparse and dense 
arrays of dots. In both cases we assume that the insulator 
separating dots is imiformly doped by donors, which do- 
nate their electrons to dots. Random positions of dots in 
the sparse case (Fig. 1) and random positions of donors 
in the dense periodic model (Fig. 5) lead to random 
charging of dots in the global ground state and to filling 
of the gap of the bare density of individual dots at the 
Fermi level. The long range Coulomb interaction creates 
a soft Coulomb gap on the background of the finite bare 
density of ground states in the vicinity of the Fermi level. 
At low enough temperature this leads to the ES variable 
range hopping conductivity, in agreement with multiple 
experimental observations. 

We concentrate on the dependence of the bare density 
of states on a number of donors per dot. Such depen- 
dence for the sparse and dense models are shown in Fig. 
(7) and Fig. (9) are qualitatively different. The sparse 
random model shows linear growth of BDOGS at weak 
doping and oscillations of the density of states at strong 
doping. On the other hand the dense model shows very 
small density of states at weak doping and no oscillations 
at stronger doping. The bare density of states determines 
the width of the Coulomb gap and the characteristic tem- 
perature of transition from ES law to Mott's law with 
increasing temperature. It determines many thermody- 
namic properties of the array as well. 

For dense arrays we have calculated characteristic tem- 
perature To of ES law in the case, when tunneling conduc- 
tance between granules is small. The challenging ques- 
tion of calculation of Tq in the case of even closer dots 
will be addressed in the next publication. Another in- 
teresting question left beyond the scope of this paper is 
the origin of ES law in isotropic and anisotropic arrays 
of strongly anisotropic (elongated) granules, for example, 
nano-tubes. 

In this paper we concentrated on ohmic hopping con- 
ductivity. However, our results can be applied to non- 
ohmic conductivity, too. It was shown^^ that at low 
temperatures a strong electric field E replaces temper- 
ature T in the exponential temperature dependence of 
the variable range hopping by the effective temperature 
Te = eE^/2kB- As a result non-ohmic ES law reads 
j = joCxp[-(£'o/£;)i/2] with Eq = 2To/eC. This depen- 
dence was observed in granular metals and nanocrystal 
thin films^'ii. 
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APPENDIX: MATRIX ELEMENTS Py FOR 
NEIGHBORING DOTS. 



Here we study the coefficients Py in Eq. (36). We know 
the diagonal term Pa ~ 1/KeffR, and for distant neigh- 
bors rij » R one can show that P^ ~ l/Kefffij- But 
the coefficients between close neighbors require further 
consideration. 

The dense array of periodic dots is equivalent to the 
lattice of identical capacitors connecting adjacent dots. 
Such a system can be studied in the language of an equiv- 
alent resistor network. Imagine a cubic lattice network 
with an identical resistance between every nearest neigh- 
bor sites. Suppose that current I goes into a lattice 
site A and travels out through the network to infinity. 
Label the potential of site A to be V ^ the potential of 
nearby dots B, C, D, • • • to be uabV, aAcV, ccadV, ■ ■ ■. 
What are the values of the uab^ olac-, c(ac, and so on? 
Here we are using the dimensionless coefficients such as 
uab = Pab/Paa, and we know all of them are between 
and 1. 

A standard numerical calculation uses the fact that 
the current going into a site equals the current going 
out of it (Kirchhoff law), which means the potential of a 
site equals the average potential of its nearest neighbors. 
Using lattices oi N x N x N sites and keeping the po- 
tential at central cite to be 1 and that of boundaries to 
be 0, one adjusts the potential values of the other sites 
so that the potential of each site is very close to the av- 
erage potential of its nearest neighbors. Using lattices 
with N up to 97 we find a = aoi = 0.34 ± 0.005, 0102 = 
0.166 ± 0.005, ao3 = 0.106 ± 0.005, an = 0.216 ± 0.005, 
ctiii = 0.17 ± 0.005 (here 111 represents the vector be- 
tween the two dots measured in the units of the lattice 
constant). 

There is also an analytic way to find the coefficients^^. 
Imagine the situation when there is charge Qq on the 
central dot while all the other dots remain neutral, one 
can use the equation similar to Kirchhoff law mentioned 
above 

r' 

where r' are 6 nearest neighbor of the site f, and (f>{r) is 
defined only on the discrete sites r = {nxR,nyR,nzR)- 
Eq. (49) can be solved using discrete Fourier transform 
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4>k~^^k- (cos kxR + cos kyR + cos k^R) = 



(50) 



One can find d)r and 



dQo f d^k 



R^ cos(fc • r) 



2kR'^ J (27r)^ 3 — cos k^R — cos kyR — cos kzR 

BZ 



(51) 



where the integration is over the Brillouin zone 
kx,ky,kz G (— 7r/i?, 7r/i?). One can use numerical in- 
tegration to find 

-Poo = = 52~' a = "01 = ^ = 0.3405, 

Qo kR^ -Poo 

ao2 = 0.1697 , Q!03 = 0.1089 , an = 0.2183 (52) 

and so on. These results agree with the results of the 
first method. 
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